Ab initio coupled-cluster study of ^^O 
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We report converged results for the ground and excited states and matter density of using re- 
alistic two-body nucleon-nucleon interactions and coupled-cluster methods and formalism developed 
in quantum chemistry. Most of the binding is obtained with the coupled-cluster singles and doubles 
approach. Additional binding due to three-body clusters (triples) is minimal. The coupled-cluster 
method with singles and doubles provides a good description of the matter density, charge radius, 
charge form factor, and excited states of a 1-particle-l-hole nature, but it cannot describe the first 
excited 0"*" state. Incorporation of triples has no effect on the latter finding. 

PACS numbers: 



One of the most important problems in nuclear physics 
is to understand how nuclear properties arise from the un- 
derlying nucleon-nucleon interactions. Recent progress 
using Monte Carlo and diagonalization |3| techniques 
produced converged results for nuclei with up to A = 12 
active particles, yielding a much-improved understand- 
ing of nuclear forces in light systems. One also must 
explore alternative methods that would not suffer from 
the exponential growth of the configuration space, en- 
abling accurate ab initio calculations for medium-size nu- 
clei. Coupled-cluster theory Q is a promising candidate 
for such developments since it provides an accurate de- 
scription of many-particle correlations at relatively low 
cost, as has been demonstrated in numerous chemistry 
applications Recently, Mihaila and Heisenberg per- 

formed coupled-cluster calculations for the binding en- 
ergy and the electron scattering form factor of ^^O using 
bare interactions In previous work 'q, we took an- 
other route and used quantum chemical coupled-cluster 
methods and the renormalized Hamiltonian to compute 
ground and excited states of "^He and ground-state ener- 
gies of ^^O in a small model space consisting of 4 major 
oscillator shells, demonstrating promising results when 
compared with exact shell-model diagonalization. 

In this Letter we report, for the first time, converged 
coupled-cluster calculations for ground- and excited- 
state energies and other properties of ^^O using mod- 
ern nucleon-nucleon interactions derived from effective- 
field theory Our ground-state calculations involv- 
ing one- and two-body components of the cluster op- 
erator are performed in up to 8 major oscillator shells 
(480 uncoupled single-particle basis states), while the 
corrections due to three-body clusters and computations 
of excited states and nuclear properties involve up to 7 
major oscillator shells (336 single-particle states). The 
significant progress in going from model calculations us- 
ing 80 single-particle states to large-scale calculations 
involving 16 correlated nucleons and almost 500 single- 



particle states has been possible thanks to the develop- 
ment of general-purpose coupled-cluster computer pro- 
grams for nuclear structure, using diagram factorization 
techniques adopted by quantum chemists. We pay par- 
ticular attention to three aspects of the calculations: (i) 
the convergence of the ground-state energy with respect 
to the size of the model space and the role of higher- 
than-two-body clusters in such studies, (ii) the ability of 
coupled-cluster methods to describe excited states, and 
(iii) the performance of coupled-cluster methods in stud- 
ies of nuclear radii, matter density, charge form factor, 
and occupation numbers. We have not yet included the 
three-nucleon interaction that should eventually be con- 
sidered However, our calculations represent a dra- 
matic step forward in nuclear many-body computations 
due to the enormous oscillator space we probe through 
application of computationally efficient coupled-cluster 
methods. They teach us about the nucleon correlations 
and the magnitude of the (missing) three-body forces. 

We use two variants of effective-field-theory-inspired 
Hamiltonians, Idaho- A and N3L0 9]. The Idaho-A po- 
tential was derived with up to chiral-order three dia- 
grams while N3L0 includes chiral-order four diagrams, 
and charge-symmetry and charge-independence break- 
ing terms. We also include the Coulomb interaction 
with the N3L0 calculations. Since very slow conver- 
gence with the number of single-particle basis states was 
obtained using bare interactions we renormalize the 
bare Hamiltonian using a no-core G:matrix approach 
which obtains a starting-energy dependence uj in the two- 
body matrix elements G{lj). We use the Bethe-Brandow- 
Petschek theorem to alleviate much of the starting- 
energy dependence (see ^3 for details). The dependence 
upon the starting energy is weak for ^^O, particularly for 
the matrix elements below the Fermi surface The 
effective Hamiltonian for coupled-cluster calculations is 
H' = t + G{uj), where t is the kinetic energy. We correct 
H' for center-of-mass contaminations using the expres- 
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sion H = H' + pc.m.Hc.m.- We choose /3c. m. such that 
the expectation value of the center-of-mass Hamiltonian 
Hem. is 0.0 MeV. We note that intrinsic excitation ener- 
gies are virtually independent of /3c. m. while the unphys- 
ical, center-of-mass contaminated states show a sharp, 
nearly linear dependence of excitation energies on /3c. m. ■ 
This allows us to separate intrinsic and center-of-mass 
contaminated states. 

Once the one- and two-body matrix elements of the 
center-of-mass-corrected effective Hamiltonian are con- 
structed, we solve the A-body problem using quantum 
chemical coupled-cluster techniques. In the ground-state 
calculations, we use the CCSD ("Coupled-Cluster Sin- 
gles and Doubles") approach to describe correla- 
tion effects due to one- and two-body clusters, and the 
CR-CCSDfT) ("Completely Renormalized CCSD(T)") 
method to correct the CCSD energies for the effects 
of three-body clusters ("Triples"). In the excited-state 
and property calculations, we use the equation-of-motion 
(EOM) CCSD method [la (equivalent to the linear re- 
sponse CCSD approach [la|)- We also correct the ener- 
gies of excited states obtained with EOMCCSD for the 
effects of triples using the CR-EOMCCSD(T) approach 
. The details of the above methods can be found else- 
where [ll[lllil|. Here, we only mention that the CCSD 
method is obtained by truncating the many-body expan- 
sion for the cluster operator T in the exponential ansatz 
exploited in coupled-cluster theory, | 'So) = exp(T) | $), 
where | ^'o) is the correlated ground-state wave func- 
tion and I <I>) is the reference determinant. The trun- 
cated cluster operator used in the CCSD calculations 
has the form T ^ Ti + Tj, where Ti = Ei,a4«Iai 

and T2 = 3 Z)jj,a6C6«I«6«j»* ^I'e the singly and dou- 
bly excited clusters and . . . (a, 6, . . .) label the single- 
particle states occupied (unoccupied) in |<I>). We deter- 
mine the singly and doubly excited cluster amplitudes 
and t^^^ by solving the nonlinear system of algebraic 
equations, {<^t\H\<^) = 0, (^^/liJl-J) = 0, where H = 
exp(-T) 7? exp(T) and |<i>°} and I'&fj') are the singly and 
doubly excited determinants, respectively, relative to |$). 
We calculate the ground-state energy i?o as (<i> | | $). 
We diagonalize the similarity-transformed Hamiltonian 
H in the relatively small space of singly and doubly ex- 
cited determinants and I'l'^j') to obtain the excited- 
state wave functions j^*^) and energies E^. The right 
eigenstates of H, where R^^'^ ^ Rq + Ri + R2 

is a sum of the relevant reference (Rq), one-body (i?i), 
and two-body {R2) components define the excited-state 
"ket" wave functions l^*^) = i?*^^) exp(T)|$), whereas 
the left eigenstates ($|iy('^^ define the "bra" wave func- 
tions {^fi\ = (^IL^^-* exp(— T). Here, each rt-body com- 
ponent of i?*^^-* with rt > is a particlc-holc excitation 
operator similar to r„, whereas L*^^^ is a hole-particle 
deexcitation operator, so that Li = ^ Ifalaa and 

-^2 = -J J2ij,ab ifj'alo'lo,baa- The right and left eigenstates 



of H form a biorthonormal set, = (5^^. If 

the only purpose of the calculation is to obtain excita- 
tion energies, the left eigenstates ($|L('^' are not needed. 
However, for properties other than energy, both right 
and left eigenstates of H are important. In particular, 
we calculate the one-body reduced density matrix p^js in 
quantum state j^*^) as follows: 

p^fi = (<i>|L('') [cxp(-T) aiap exp(r)] |$) . (1) 

In the CCSD ground-state (/i = 0) case, we have 
T ^ Ti + T2, = 1, and L(o) = 1 -h Ai -h A2, 

where the one- and two-body deexcitation operators Ai 
and A2 are determined by solving the CCSD left eigen- 
value problem, obtained by right-projecting the equation 
I {1 + A)H = Eo{^\{l + A), with Eq representing 
the CCSD energy and A = Ai -I- A2 on the singly and 
doubly excited determinants. Thus far, we have focused 
on the CCSD and EOMCCSD methods which use inex- 
pensive computational steps that scale as n^n^, where 
n-o (riu) is the number of occupied (unoccupied) single- 
particle states. While the full inclusion of triply excited 
clusters is possible, the resulting methods are expensive 
and scale as n^n,^. Thus, we estimate the effects of T3 
and i?3 on ground- and excited-state energies by adding 
the corrections to the CCSD/EOMCCSD energies, which 
only require n^n^ noniterative steps. These corrections, 
due to T3 and R3, define the CR-CCSD(T) and CR- 
EOMCCSD(T) approaches |i,[l3. In this study, we use 
variant "c" of the CR-CCSD(T) and CR-EOMCCSD(T) 
approaches 0. 

We turn to a discussion of our ^^O results. We 
choose the oscillator energy Tiuj for our basis states 
to minimize the CCSD energy. For the N — 7 
and = 8 oscillator shell runs, Tiuj = 11 MeV, 
and the results are nearly independent of hto |lOl |. 
Shown in Fig. □ are our CCSD/EOMCCSD and CR- 
CCSD(T)/CR-EOMCCSD(T) ground- and excited-state 
energies as a function of N. The symbols in Fig. ^ rep- 
resent our calculations while the lines represent a fit of 
the form E{N) — E^o + aexp {—bN), where the extrap- 
olated energy i?oo and a and b are parameters for the 
fit. We also show in Fig. ^ our calculations for the first 
excited 3~ state and the position of the lowest calculated 
0"^ excited state. We now discuss these results. 

Triples correction to the CCSD ground-state energy. 
The small model space calculation 7] indicated that the 
triples corrections to the ground-state CCSD energies are 
small. We extended these calculations from 4 to 8 major 
oscillator shells for CCSD calculations and to 7 major 
oscillator shells for CR-CCSD(T) calculations, as shown 
in Fig. n We find that the extrapolated CCSD energy is 
-119.4 MeV for Idaho-A. For the = 7 Idaho-A calcula- 
tion, the difference between the CCSD and CR-CCSD(T) 
result is 0.6 MeV, while the extrapolated values differ 
by only 1.1 MeV; our extrapolated CR-CCSD(T) en- 
ergy is —120.5 MeV. The Coulomb interaction adds to 
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the binding 11.2 MeV, so that our estimated Idaho-A 
ground state energy is —109.3 MeV (compared to an ex- 
perimental value of -128 MeV). Our N ^ 7 {N = 8) 
N3L0 CCSD and iV = 7 CR-CCSD(T) energies, which 
include the Coulomb interaction, are —112.4 (—111.2) 
and —112.8 MeV, respectively. Thus, the two-body in- 
teractions underbind "'^^O by approximately 1 MeV per 
particle, pointing to the need for three-body forces. For 
the Idaho-A and N3L0 interactions and the ^^O nucleus, 
we conclude that connected T3 clusters are indeed small, 
contributing less than 1% to the ground-state energy. 
This is an important finding, since it implies that es- 
sentially all correlations in a closed-shell nucleus result- 
ing from two-body nucleon-nucleon interactions can be 
captured by the relatively inexpensive CCSD approach. 
Another important finding is a rapid convergence of the 
CCSD and CR-CCSD(T) energies with the number of os- 
cillator shells owing to the use of the renormalized form 
of the Hamiltonian. For example, the difference between 
the iV = 8 and iV = 7 CCSD/Idaho-A energies is 0.5 
MeV (see Fig.Hl). 

Calculations of the first excited 3~ state. The first 
excited 3~ state in ^^O is thought to be principally a 
one-particle-one-hole (Ip-lh) state The experience 
of quantum chemistry is that the EOMCCSD and CR- 
EOMCCSD(T) methods describe such states well, pro- 
vided that the three-body interactions in the Hamilto- 
nian can be ignored. The largest i?i amplitudes obtained 
in the EOMCCSD calculations indicate that the domi- 
nant Ip-lh excitations are from the Opi/2 orbital to the 
Od5/2 orbital. The 2p-2h excitations in the EOMCCSD 
wave function, defined as R2 + RiTi + Ro{T2 + Tl/2) 
(i?o = in this case), are much smaller than the Ri am- 
plitudes, and the CR-EOMCCSD(T) calculation hardly 
changes the total energy of the state, which indicates 
that this state has indeed a Ip-lh nature. Our ex- 
trapolated Idaho-A results indicate that the 3~ state 
lies at -108.2 and -108.4 MeV in the EOMCCSD and 
CR-EOMCCSD(T) calculations, respectively The CR- 
EOMCCSD(T) method yields an excitation energy of 
12.0 MeV for this state which experimentally lies at 
6.12 MeV. N3L0 yields similar results. Based on the 
Ip-lh structure of the state, we conclude that Idaho-A 
and N3L0 do not yield an excitation energy for the 3~ 
state which is commensurate with experiment. These re- 
sults agree with recent no-core shell-model calculations 
with similar two-body Hamiltonians The 3~ state 

is expected to be built on Ip-lh excitations which de- 
pend on the single-particle splittings. These splittings 
will be affected by three-body forces not included in our 
Hamiltonian, thus affecting the energy of the 3^ state. 
Whether other mechanisms than three-body forces can 
provide an additional binding of 6 MeV needs further re- 
search. Our results are converged at the coupled-cluster 
level employing the Idaho-A and N3L0 two-body in- 
teractions, so it is likely that the discrepancy between 



theory and experiment resides in the Hamiltonian, not 
in the correlation effects which EOMCCSD and CR- 
EOMCCSD(T) describe very well if three-body forces 
play no role and if the state has a Ip-lh nature. 

Calculation of the first excited 0+ state. This state 
(experimentally at 6.05 MeV), believed to have a Ap-Ah 
character, cannot be described by EOMCCSD or CR- 
EOMCCSD(T). This is confirmed by our calculations as 
we see large differences between the EOMCCSD or CR- 
EOMCCSD(T) results and experiment (see Fig.[lJ. One 
would need to include Ap-Ah operators (T4 and R4) to 
improve coupled-cluster results. 
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FIG. 1: The energies of the ground-state (g.s.) and first- 
excited 3~ and 0"'' states as functions of the number of oscil- 
lator shells A'' obtained with coupled-cluster methods and the 
Idaho-A interaction. 

Although we concentrated on the lowest energy J = 3~ 
and J = 0''" excited states and the role of three-body clus- 
ters on these, we also performed preliminary calculations 
for other negative parity states. The quartet of negative 
parity states starting with the J = 3~ state, and includ- 
ing the J = 1~, 2^ and 0~ states, are all believed to have 
a similar Ip-lh character [l3|. The EOMCCSD calcula- 
tion with 5 major oscillator shells and Idaho-A confirms 
the existence of this quartet, giving excitation energies 
of 13.57, 15.37, 17.07, and 17.15 MeV for the J = 3", 
1~, 2~, and 0^ states, respectively. While these states 
are all a few MeV above the experimental values, their 
ordering predicted by EOMCCSD is correct. 

Calculation of the one-body density. We use Eq. 
where /i = 0, to calculate the ground-state density for 
^^O. We show the resulting radial density, p{r), in Fig.|2 
The root-mean-square (rms) radius is found through an 
integration r'^.^^ — J r'^p{r)dr/ J r'^p{r)dr. To obtain a 
charge radius, we correct this value for the finite size of 
the nucleons, which experimentally are = 0.743 fm^ 
and = 0.115 fm^, and for the Os center-of-mass mo- 



62.2071 

Ahuj 



fm^. Our 



tion, for which we use (^'o | R | ^"0) 
rms charge radii for ^^O for 5, 6, and 7 oscillator shells 
are 2.45 fm, 2.50 fm, and 2.51 fm, respectively when the 
Idaho-A potential is used (N3L0 gives similar values). 
The experimental charge radius is 2.73±0.025 fm. We 
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FIG. 2: Top panel: The charge form factor computed from 
the CCSD density matrix. Bottom panel: the matter density 
in ^^0. The results obtained with the Idaho-A interaction. 



lowest J = 3~ state, which is, quite likely, due to an in- 
adequate description of the relevant nuclear forces by the 
Hamiltonian. We were unable to accurately describe the 
lowest J = 0+ excited state due to connected Ap-Ah cor- 
relations missing in coupled-cluster approximations em- 
ployed in this study. The CCSD method provides reason- 
able results for the nuclear matter density, charge radius, 
and charge form factor. The use of the renormalized 
Hamiltonian guarantees fast convergence of the results 
with the number of oscillator shells. All of this makes 
low-cost coupled-cluster methods a promising alternative 
to traditional shell-model diagonalization techniques. 
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also calculated the occupation probability for the natu- 
ral orbitals. Experimental data from quasi-elastic proton 
knockout HI yields 2.17 ± 0.12% for the Odg /2 occupa- 
tion and 1.78±0.36% for the lsi/2 occupation. We obtain 
3.2% and 2.3% respectively, using Idaho-A in the N = 7 
model space. For N3L0 in the N = 7 model space, we 
obtain 3.8% and 2.6%, respectively. For the calculation 
of the nuclear charge form factor, we follow [23|. In this 
approach, the form factor includes contributions from the 
two-body reduced density matrix due to centcr-of-mass 
corrections. We computed the one-body density con- 
tributions within the framework of CCSD theory using 
Eq. 1^. The contributions of the two-body density ma- 
trix were computed within the shell-model like descrip- 
tion as = (^0 I a'lapasaj \ '^o)/ (*o|*o), where we 
approximated | *o) by (1 -I- Ci -I- C2) | with Ci = Ti 
and C2 = T2 -I- ^T^ defining the Ip-lh and 2p-2h com- 
ponents of the CCSD wave function. The upper part of 
Fig. 121 shows the charge form factor for different model 
spaces. The 5-shell and 6-shell results include the center- 
of-mass corrections and exhibit a second zero. Compared 
to the experimental value (the arrow in Fig. [21), the first 
zero of the form- factor is reasonable, although slightly 
too large; this is consistent with an underestimated value 
of the theoretical charge radius. 

In summary, the ^^0 ground state is converged with re- 
spect to the model space size and is accurately described 
within the basic CCSD approximation, with three-body 
clusters contributing less than 1% of the binding energy. 
We attribute the 1 MeV per particle difference between 
the coupled-cluster and experimental binding energies to 
three-body forces. We obtained a correct description of 
the quartet of low- lying negative parity Ip-lh excited 
states, although there is a 6-MeV difference between the 
converged coupled-cluster results and experiment for the 
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